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Abstract 

We analyze numerically the Balitsky-Kovchegov equation with the full impact pa- 
rameter dependence b. We show that due to the particular ^-dependence of the initial 
condition the amplitude decreases for large dipole sizes r. Thus the region of satura- 
tion has a finite extension in the dipole size r, and its width increases with rapidity. 
We calculate the ^-dependent saturation scale and discuss limitations on geometric 
scaling. We also demonstrate the instant emergence of the power-like tail in impact 
parameter, which is due to the long range contributions. Thus the resulting cross 
section violates Froissart bound despite the presence of a nonlinear term responsible 
for saturation. 



1 Introduction 



The high energy hmit of QCD is one of the most intriguing aspects of hadronic physics. 
With the advent of new generation accelerators hke HERA, Tevatron, RHIC and in a near 
future LHC, the basic problems of strong interactions are experimentally studied and con- 
fronted with theoretical predictions of high energy QCD. The important discovery of the 
rise of the proton structure functions at HERA 1 at small values of the Bjorken variable 
X (which is equivalent to the high energy limit) is an example of such a confrontation. 
This rise was predicted by high energy QCD [2 and is related to the increase of the gluon 
density. Ultimately the rise has to be damped by the presence of the saturation effects 
which enter via nonlinear modification to the QCD evolution as first proposed in the pio- 
neering work . Thus QCD at high energy is the theory of high density systems of colored 
particles. In such systems hard scales appear which allow to apply perturbative techniques 
although effective interactions in the dense partonic system are of non-perturbative origin. 
The interplay between hard and soft (perturbative and non-perturbative) aspects of QCD, 
which also touches in the operational way the issue of confinement, is the most exciting 
element of high energy QCD. 

The effective theory which describes high energy scattering in QCD is Color Glass 
Condensate . The basic equations of this theory |5] are equivalent to the hierarchy of 
equations derived by Balitsky [Jj and later on reformulated in a compact form by Weigert 
[S]. These equations contain the BFKL evolution and also the triple Pomeron vertex 
in]. In this work we present the results of numerical studies of these basic equations. To 
be precise we study the simplified version of the hierarchy of Balitsky's equations which 
reduces to one equation in the limit of large number of colors. This is the equation 
obtained by Kovchegov [lUj in the dipole approach IllJ to high energy scattering in QCD. 
With this equation the deep inelastic lepton-nucleus scattering can be described and also 
information on small-x hadronic wave function be obtained. 

We study this equation in the full form, including the impact parameter dependence b. 
Previous analytical ^21 El and numerical [m El E] studies were done under a simplified 
assumption of infinitely large and uniform nucleus, i.e. neglecting impact parameter b. 
Recently, an approximate solution to this equation in semiclassical approach was consid- 
ered ^7]. As we will show, the solutions of the full form of the Balitsky-Kovchegov (BK) 
equation possess important new features in comparison to the uniform case, e.g. restricted 
scaling properties with impact parameter dependent saturation scale. The detailed anal- 
ysis of the 6— dependence of the solution allows to study the high energy behaviour of the 
7*A^ cross section. We show that the Froissart bound JHl is violated due to the long range 
contributions in the kernel jl9j which brings the issue of the lack of confinement effects in 
the BK equation. The problem of the Froissart bound was also extensively discussed in 
the same context in jl9l l2Uj based on analytical consideretations. 

The paper is organised in the following way. In the next section we briefly present the 
BK equation and its symmetries. In Sec. 3 we describe the numerical methods of flnding 
the solution and discuss the initial condition. In Sec. 4 we present the resulting amplitude 

as a function of dipole size r and extract the b - dependent saturation scale. In Sec. 5 
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we discuss the form of the impact parameter profile which emerges in the evolution, and 
in particular we concentrate on the emergence of the power tails in b in the amplitude. We 
also present the estimate of the cross section of the black disc radius and its dependence 
on the rapidity. Finally, in Sec. 6 we state our conclusions. 



2 The Balitsky-Kovchegov equation 

The deep inelastic scattering of a lepton on a nucleus at high energy in the dipole picture 
im I21j is viewed in the nucleus rest frame as the splitting of an exchanged virtual photon 
into a qq dipole and the subsequent interaction of the dipole with the nucleus. The latter 
process is described by the dipole-nucleus scattering amplitude A^(x, y), where x,y are 
two-dimensional vectors of the transverse position of the dipole ends. Alternatively, one 
can introduce the dipole vector r = x — y, and the impact parameter b = (x + y)/2. Thus 
in general, the amplitude depends on the four transverse degrees of freedom and rapidity, 
Y = ln(l/x), playing the role of the evolution parameter 

iV(x,y,y) ^ N^yiY). (1) 

From now on, for shortness of the notation, we assume the Y— dependence implicit. 

In the leading logarithmic approximation, the dipole scattering amplitude obeys a 
nonlinear evolution equation derived by Balitsky and Kovchegov El 

dN^y _ fcFz (x-y)2 



dV 



/ ^ ( ^^2r^^ ^2 {^xz + iVyz - iVxy " iVxz iVyz} , (2) 

J 27r (x - z)^(y - z)^ 



where Og = agNc/TT. In addition, one has to specify an initial condition at y = Yq: 
A'xy = A^'^(r,b). The amplitude N{x,y) in ^ is given by the following correlator 

iV(x,y) = i-Tv(l-[/t(x)C/(y)) , (3) 

where the trace is done in the colour space, and the eikonal factor U is defined as the path 
ordered exponential with the SU (N) gauge fields (in the gauge A~ = 0) 



The averaging (...) in ^ is performed over an ensemble of classical gauge fields. In general, 
an infinite hierarchy of equations is found for correlators of the U factors [7]. In the large 
Nc limit, however, the closed form (j2} can be found for the two point amplitude A'xy 
jlOj . The linear part of ((21) corresponds to the dipole version 11 of the BFKL equation 
|22j at nonzero impact parameter and its solution has been studied in the Monte Carlo 
simulation of onium-onium scattering |23| . The additional quadratic term emerges due 
to the summation of the multiple interactions of the dipoles in the quark-antiquark wave 
function with the nucleus. 



3 



Figure 1: The triangle geometry in Eq. The points x, y,z are the dipole ends in the 
transverse space, and the vectors bxz,bxy,bzy are the impact parameters of the three 
dipoles. 



The r.h.s of Eq. ^ has a nice geometrical interpretation: the parent dipole (x, y) splits 
into two new dipoles (x, z) and (y, z), and the summation is taken over all new dipoles. On 
the other hand, the nonlinear term describes the recombination of the two dipoles (x, z) 
and (y, z) into one (x, y). The three dipoles form a triangle, shown in Fig.^ The first part 
of the integral kernel in ((2)) only depends on the triangle sides: |x — y|, |x — z| and |y — z|. 
A nontrivial dependence on the position of the triangle in the plane, i.e. the dependence 
on the impact parameter vectors bxy, bzx, by^, is introduced through the arguments of the 
amplitudes A^. In particular, it is interesting to study how the 5-dependence introduced 
by an initial condition A^'^(r,b) propagates with increasing Y. Let us also note that the 
singularities at z = x, y in Eq. Q are integrable provided 

lim Axy ~ [x - y|^ e > . (5) 

x^y 

Thus the dipole which shrinks to a point does not scatter. 



2.1 Symmetries of the BK equation 

The BK equation has a rich symmetry structure. Introducing the complex number 
notation for the transverse vectors, e.g. x = xi + i X2 and x = xi — i X2 for x = (xi,X2), 
it can be easily shown that the measure in Eq. @, 



as (x - y)2 2 



27r (x-z)2(y-z)2 
is invariant under the Mobius transformation^ 



d^z , (6) 



^Provided as is kept constant. 



ax + b _ a X + b 
cx + d cx + d 
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where the parameters a,b,c,d € C and ad — be ^ 0. Identical transformations are also 
applied to y (y) and z (z). Thus the BK equation is covariant with respect to the Mobius 
transformation. In particular, the following elementary transformations from ^ are rel- 
evant for our discussion 

- global two-dimensional translations by vectors b: x — > x + b , 

- global two-dimensional rotations by angles <j): x ^ O(0)x, 

- scale transformations with a real, positive parameter A: x — > Ax 

- inversion (in complex notation) : x — > 1/x . 

Notice, that if an initial condition A^^ is invariant under any of the discussed trans- 
formations, the solution of the BK equation N^y{Y) preserves the corresponding sym- 
metry. In particular, if an initial condition is invariant under translations and rotations, 
A^(x,y) = iV''(|x — y|), the solution at any rapidity Y has the same property, i.e. it only 
depends on the dipole size r = |x — y| but not on the impact parameter b = (x + y)/2. 
The problem of finding solution simplifies enormously in this case since only one degree 
of freedom is relevant, namely the dipole size r. 

Physically, this approximation (called local approximation) corresponds to an infinitely 
large and uniform nucleus. Previous analytical jl21 113j and numerical jl4l 1151 116j studies 
of the BK equation were based upon this assumption. In this approximation, the solution 
shows saturation, N{r) — > 1, with the characteristic scale Qs{Y). For dipoles smaller than 
the inverse of the saturation scale, r < l/Qs{Y), the solution is governed mainly by the 
linear term of Eq. ^ and shows the exponential rise in rapidity, ~ exp(a;jpy) where 
wjp = 4 In 2 is the intercept of the BFKL kernel. On the other hand, in the region where 
dipoles are large, r > l/Qs(Y), the nonlinear term slows down the rise and eventually the 
amplitude saturates to 1. The saturation scale QsiY) depends on the rapidity in the 
following way m Ha US] 

Qs(Y) = Qoexp{XasY), (8) 

where the coefficient^ A ~ 2. The solution in the local approximation also exhibits a 
property of the geometric scaling namely for r > l/Qs{Y) 

Nir,Y) = NirQsiY)) , (9) 

which means that the amplitude in the saturated region only depends on one combined 
variable rQs{Y) instead of r and Y separately^. Thus the diffusion into infrared, typical 
for the linear BFKL equation, is damped by the emergence of the saturation scale Qs(Y) 

m 

^To be precise in Ref. |12) the coefficient was found to be tlie same as tlie Pomeron intercept A = 
ujp/as = 41n2, but this value was not confirmed by subsequent analytical and numerical studies. 
^This is also a feature of the saturation model |27|. 
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Figure 2: Dipole degrees of freedom: the dipole vector r and the impact parameter b. 

Since in our analysis we want to study the impact parameter dependence, we obviously 
have to abandon the assumption about the translational invariance. However, in order to 
simplify the problem, we adopt more physical assumption that our nucleus is cylindrically 
symmetric, i.e. A^(x, y) is invariant under the global rotation. It means that in the 
parameterisation of the dipole position, see Fig. [21 

(x,y) = ir,b,e,cP), (10) 

we drop the dependence on the azimuthal angle (p. Note that we keep the dependence on 
6 which is the angle between the vectors r and b 

r • b 

cosO = (11) 
r 

The assumption about cylindrical symmetry reduces the number of parameters in the 
amplitude to four: three degrees of freedom for a dipole and the evolution variable Y. 

3 The numerical method of finding solution 

The numerical method for the solution of the BK equation is very similar to the one used 
for the solutions of the linear equations, see |55j. We discretise the amplitude N{r, b, cos 6) 
in all 3 variables (Inr, ln6, cos^). A simple linear interpolation has been chosen which is 
the fastest in this case. To find the solution as a function of rapidity we take in the first 
step , the initial condition at y = lo> and evaluate Eq. @ with the step Ay. This 
gives the first approximation for the solution at y = Ay 

My^ = + AY [ {iVO, + iVO, - iVOy - AtO^AtO, , } (12) 

where we symbolically denoted the integration over z with the measure ©. If the relative 
difference |(A^4y - A^xy)/^xyl 

< e (some accuracy) at each point of the grid, we finish our 
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procedure, otherwise solution (|12() is used to find the second approximation from: 




The r.h.s. of Eq. (|13() was found after integrating over Y from to AF, assuming a linear 
interpolation in Y between A^*^ and N^^\l\Y). This is justified for small enough value of 
Ay. We iterate in this way, using Eq. (|13|) . until the desired accuracy is achieved. 

Usually only couple of iterations are needed to find the right answer, in fact we fix the 
maximal number of iterations to 5. In order to get satisfactory results one has to work 
with a grid which is at least (100,. x 100^ x 20c). Each step in rapidity produces about 
1.5 MB of output data and takes about 300 min. of CPU time when run on PC machine 
with 2.5 GHz processor and 2.0 GB RAM memory. 

One of the important issues while studying the BK equation is the choice of the initial 
condition at rapidity y = 0. As mentioned above, because of the cylindrical symmetry 
our initial amplitude A^'' should only depend on the three variables (r, 6, 6). Since we have 
nearly no information on the angle between r and b, we do not assume any dependence 
on it in the initial conditions. We shall see that the nontrivial dependence is nevertheless 
generated through the evolution. As far as the r and h dependence is concerned we have 
chosen the initial distribution in the Glauber-Mueller form |29l l30j 

iV°(r, 6) = 1 - exp{-r25(6)} , (14) 

with S{h) being a steeply falling profile in 6, e.g. S{h) ~ exp (-6^). The Glauber-Mueller 
formula resums the multiple scatterings of a single dipole on a nuclear target and it has 
been advocated to be a natural choice for the starting distribution of the BK equation 
[m] . Let us note that the above distribution has a property that for any fixed value of 
impact parameter it saturates to 1 for sufficiently large values of the dipole size r. More 
detailed analysis with other forms of the initial conditions will be presented elsewhere |31) . 

4 The dependence of the solution on a dipole size r 

We have performed the numerical evolution of the BK equation starting from formula H14|l 
as the initial condition with the following impact parameter profile 

5(6) = 10 exp{-6V2} . (15) 

Throughout this work we keep the coupling constant in @ fixed, = 0.2. The running 
of the coupling, although physically more justified, would be an additional complication 
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Figure 3: The amplitude N(r, b, 0; Y) as a function of dipole size r for indicated values of 
rapidity Y and two values of the fixed impact parameter: b = 0.2 (left) and b = 5 (right). 
The orientation of the dipole cos 9 = 0. The red-dashed line is the input distribution H14|) 
with profile (fT3|l. 

to the solution which we would like to avoid at this stage^. Nevertheless, one should 
stress that a running coupling is an important NLL effect [H3], which should be taken into 
account in future investigations. We are going to study now how such initial condition 
(|14j) is modified by the evolution in rapidity. We concentrate first on the dependence of 
the solution on a dipole size r. 

The results of our numerical analysis are shown in Fig. |21for two values of the impact 
fixed impact parameter: b = 0.2 and 6 = 5. We have also fixed the orientation of the 
dipole such that r _L b that is cos^ = 0, see Eq. As we see, for the values of r 

which are small (here r <^ 1) the amplitude rises with increasing values of rapidity Y. In 
this region, the linear part of the BK equation dominates. With increasing r, for fixed 
rapidity, the amplitude finally reaches the saturation value A'sat = Ij were the nonlinearity 
of the BK equation is crucial. This is mostly visible for small value of impact parameter. 
Fig. 13 (left). For a larger value of impact parameter, 6 = 5, in the plot to the right, the 
saturation has been only reached for the highest value of rapidity shown Y = 11. What 
is interesting is the fact that the amplitude has a maximum for the dipole size which is 
twice its impact parameter r = 26, that is r = 10 in this case. 

At large r, the behaviour of the amplitude is quite different from the case with the 
translational invariance (no 6— dependence). The amplitude decreases with increasing r 
even if the initial condition is saturated {N^ = 1) in this region. The reason for this 

^In the local approximation, running coupling changes rapidity dependence of the saturation scale: 
QT'iY) ~ exp(c/F), see [11101103 and also 0. 
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behaviour can be easily understood analysing the solution after the first iteration ()12|) 



nS = K + £/ {K + - < - N^y.} , (16) 

with the Glauber-Mueller initial condition 

ArO(r,6) = I - exp{-cr'^ e^p{-b'^d)} . (17) 

Saturation means that for a large enough r = |x — y| (and fixed b) N^y — > 1, i.e. a large 
dipole is totally absorbed. This is indeed the case for the initial condition H17|) . After the 
first evolution step, however, A^y becomes less than 1 for large r due to the large and 
negative integrand in the curly brackets in H16|). which leads to the effect shown in Fig. |21 
at large r. 

Indeed, if N^y ~ 1 in Eq. the integrand is always negative or equal zero 

+ N^.-l- iV°,iV°, = -(1 - iV°y)(l - A^°J < . (18) 

In Fig.Qwe show an example of typical configurations of the (x, z) and (y, z) dipoles which 
give dominant contribution to the integral in Eq. ()16j) . For such a contribution rxz ~ ~ 
2bxz ~ 2byz ~ r (for bxy ~ 0), and for the initial condition (fTTj) the amplitudes N^^ k, 
Ny^ ~ for the sufficiently large value of r. Thus the integral picks up significant negative 
contribution from this configuration. This should be contrasted to the translationally 
invariant case with no 6— dependence in the initial condition, e.g. N^{r) = 1 — exp{— r^}. 
For such initial condition the configuration from Fig. ^ gives 

^xz ~ A'yz - 1, and the 

expression in Eq. (|18j) is negligible. One can prove that any other configuration, for 
example when |x — z| <^ |x — y| or |x — z| ^ |x — y|, leads to the same result. Thus, in 
the case of infinitely large and uniform nucleus, the amplitude N stays always saturated 
for large dipole sizes. 

In the next evolution step, plays the role of the initial condition. S ince Axy is no 
longer equal to 1 for large dipole sizes, the analysis becomes more complicated. At some 
rapidity, the discussed integral becomes positive leading to the effect observed in Fig. |21 
where N starts to grow again, reaching unity at large r for large values of rapidity. 

In more physical terms the fall-off of the amplitude at large r can be explained by 
realizing that the end points x and y of a sufficiently large dipole are in the region where 
there is no gauge field. In this case ?7(x) = U{y) = 1, and the correlator Q vanishes. It 
simply means that the dipole is so large that it misses the localised target given by S{b). 
On the other hand, the non-vanishing value of the amplitude for very large values of r in 
the previous studies in the local approximation ^1 El El El Ej was a consequence of 
the fact that the field was uniform and present everywhere, thus it did not matter how 
large the dipole was, it always scattered. The appearance of the gauge field, signalled by 
the increase of N with rising rapidity, corresponds to expansion of a black (or grey) region 
of a nucleus, see section 15.31 
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4.1 Saturation scale and geometric scaling 

The saturation scale Qs is defined from the condition 

^N{r = l/Q,Ae,Y)), = K, (19) 

where the constant k is of the order of unity {k = 1/2 in our numerical analysis) and 
we have taken an average over the orientation of the dipole with respect to the vector 
of impact parameter. Notice that after solving the above equation, the saturation scale 
depends on impact parameter in addition to rapidity^. For the form of N , shown in Fig.|51 
there are two solutions of Eq. H19|) for sufficiently large rapidities, corresponding to small 
and large dipole sizes. Thus saturation (defined here as, N > k) occurs over a finite region 
of dipole sizes 

l/Qs{b,Y)<r<RH{h,Y). (20) 

The scale Qs{h,Y) is the 6-dependent saturation momentum introduced by Mueller 
[30] . The 6— independent saturation scale found in the previous analyses [T3 l IT U ITb I ITB] 
can be regarded as an average over all area of interaction QsiY) = {QgiY^ b))b- '^^^ second 
scale Rnib^Y) appears because we have introduced exponential impact parameter profile, 
as has been explained in detail in the previous section, and refiects the boundary of the 
nucleus. In Fig. |l]we show the form of the impact parameter profile of the saturation scale 
Qs{b,Y) for two different rapidities Y = 5 and Y = 11. The saturation scale at y = 0, 
shown by the dashed line for comparison, is computed from the initial condition 1)14(1 with 
(fT5|) 

Q^(6,0) = Q2exp(-6V2), (21) 

where we have fixed the normalisation Qo to match the Qs{b, Y) at small values of 6 = 0.1. 
Thus the dashed lines in Fig. 4, show the profile of Qs{b,Y) which we would get if the 
exponential impact parameter dependence set in the initial conditions would be preserved 
through the evolution. We clearly see that for higher impact parameters, the exponential 
fall-off set by initial condition is replaced by the power-like tail Qs{b,Y = 11) ~ l/b"^ with 
7 ~ 1.6 - 2.0 for b>7. 

As it is evident from Fig. 13 both the saturation scale Qs(b,Y) and RnibjY) rise with 
increasing rapidity, thus leading to broadening of the region in r for which saturation 
occurs. The energy dependence of Qs{b,Y) has been estimated to be proportional to 
exp{XsasY) with ^ 1.7—2.0 at the highest values of rapidity Y = 9—11. We also checked 
that this value of As is not very sensitive to impact parameter. This would mean that (at 
high rapidities) the saturation scale has a factorised form Qs{Y, b) = ex.p{XsasY) g{b). We 
expect exponential dependence on rapidity of RH{b,Y) ~ exp{XHasY) too, and we have 
found Xh — 2.9. 

Since there are two scales in the problem, Qg and Rh, geometric scaling of the form @ 
is not present anymore for arbitrarily large dipoles. However, when r <^ RH(b,Y) there 
is still an approximate geometric scaling in the combined variable rQs{b,Y) for fixed b. 

^Of course, Qs depends also on k but this dependence is irrelevant for our discussion. 
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Figure 4: The dependence of the saturation scale Qs{b,Y) on impact parameter for two 
fixed values of rapidity. For the comparison the saturation scale at initial rapidity (|21|) is 
shown by the dashed line, normalised to the value of Qs{b, 1") at 6 = 0.1. 

This is shown in Fig. [51 obtained from Fig. after rescaling r by exp(1.7asy), i.e. by the 
rapidity dependence of the saturation scale Qs(b, Y), extracted in the region 7 <Y < 11. 
As discussed, there is no geometric scaling for large dipole sizes. The problem of scaling in 
the 6— dependent case has already been addressed in [20] and in recent phenomenological 
studies j3S|) and still needs deeper analysis. 

Let us finally note that for our solution there is a region in b in which the saturation 
scale does not exist at all. It is evident from the right hand side plot in Fig. 01 where up to 
the rapidities y ~ 6 the amplitude N < k for all values of the dipole sizes r, and Eq. ((T^ 
does not have a real solution. 

5 Impact parameter dependence of the solution 

We will discuss now the impact parameter dependence of the solution to the BK equation 
with the Glauber-Mueller input distribution (|14I15|) . 

In Fig. ini we show the amplitude N{r, b, 0; Y) as a function of impact parameter b for 
fixed and small dipole size r = 0.1 and various values of rapidity. We have also fixed the 
orientation of the dipole such that cos ^ = 0. The first striking feature of the solution is 
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r exp(1 .7 ttg Y) 

Figure 5: Geometric scaling for b = 0.2 obtained after rescaling r by the rapidity depen- 
dence of the saturation scale Qs{b,Y). Dotted, dashed and solid lines are for Y = 11,9,7 
respectively. The orientation of the dipole has been fixed such that cos 6* = 0. 

the fact that the steeply falling exponential dependence in 6, given by the initial profile 
S{b), is washed out by the evolution and instead clear power behaviour is generated for 
large values of b. Initially, at lowest rapidity Ay = 0.1 the amplitude N(b) ~ 1/6^'^, 
whereas later on it has a milder dependence: N{b) ^ 1/b"' with 7 between 2 and 3 for 
y > 5. The growth of the amplitude as a function of rapidity at large impact parameters 
is exponential, N{Y) ~ exp(a;y) with uj ~ 2.7 for b = 10. This strong dependence 
is clearly governed by the linear BFKL part of the equation since the amplitude is small 
enough for the nonlinear term to be safely neglected and it is in a very good agreement 
with the hard Pomeron intercept ujjp = 41n2as = 2.77 ag- 

On the other hand, at small values of impact parameters, b < 1, where the amplitude 
is large, we clearly observe that the growth of amplitude is strongly damped due to the 
nonlinear term. This is the region of b where the saturation sets in first. 

5.1 The origin of the power-like tail 

It has been argued by Kovner and Wiedeman in ^H] that the power behaviour in impact 
parameter N ~ l/b'^ originates from the configurations with very large dipoles. Following 
jl9j let us take small dipole with size r = |x — y| which is located far away from the target, 
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Figure 6: The amplitude N(r, b,Y) as a function of impact parameter 6 for different values 
of rapidity Y. The dipole size and orientation are fixed, r = 0.1 and cos 9 = 0. The dashed 
line is the input distribution ((TUl with the profile lfTB|l. 



at large impact parameter 6, in the area where the colour field is very weak. In that case 
one has f/(x) ~ U (y) ~ 1 and thus A^(x, y) ~ 0. 

The non-vanishing contribution to the r.h.s of equation comes from configurations 
of large dipoles with one end-point situated at x, and the other at z, close to the center of 
the target where the field is strong. The phase in eikonal (HJ oscillates strongly and thus 
{U{z)) ~ 0. Therefore iV(x, z) ~ N{y,z) ~ 1, and this configuration gives 

/ "^'^ (x - z)^(y - z)2 (^^^ + - - ^^^^y^^ ^ ^ / ^ F 

where 7ri?g(y) is the area of strong field in the target over which we integrate. Strictly 
speaking the statement that A^(x, z) = 1 for the configuration considered above is valid 
only at very high rapidities. Due to the initial conditions (|14() at intermediate rapidities 
there will be always such b, large enough, for which these configurations will have A^(x, z) < 
1, however always N{x,z) ^ A^(x, y). 

In order to check this statement numerically we perform only one iteration of the BK 
equation with very small step in rapidity Ay = 0.1, Eq. (|12j) . and divide the integration 
region into two parts: |z — b| < tq and |z — b| > tq with the cutoff tq = 1 for the choice 
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Figure 7: The r.h.s. of the BK equation ((2j computed for the initial condition (|14|15|) 
(solid line) and the contributions from the short and long dipoles (dashed lines), see ()23() . 
plotted as a function of b for fixed dipole size r = 0.1 and orientation cos 6 = 0. 



of small dipole r = 0.1. To be precise we evaluate two integrals: 



short 



long 



e(ro-|z-b|) + / G(|z-bl-ro) 



d?z (x — y)^ 
(x-z)2(y-z) 



' V xz ' v: 



xy 



xz yz 



(23) 

with iV° being our initial condition. The results are presented in Fig. [71 where the short 
range contribution dominates at small values of impact parameters and the long range one 
appears to be responsible for the power behaviour in b at large values. The point at which 
the long range contribution starts to dominate over the short range one is determined by 
the form of the initial condition. 



5.2 Angular dependence for large dipole sizes 

It is interesting to study also the angular dependence of the solution A^. For sufficiently 
small dipole sizes r <C Rh, the angular dependence is negligible. However, when r ~ Rh 
the difference in the amplitude due to the dipole orientation is quite substantial, especially 
at large impact parameters 6 ~ r. 
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Figure 8: The amplitude N(r,b,9;Y) as a function of impact parameter b for two rapidities: 
Y = 2 (left) and 1" = 8 (right). The dipole size is fixed to r = 20. The solid blue line 
corresponds to the case r || b, and for the dashed blue line r _L b. The red-dotted line is 
the input distribution ((Tl)) with ((T3|l . 

In Fig. IHl we plot N(r, b, 9; Y) as a function of b for a large dipole with r = 20. Two 
rapidities were considered. In both cases we have compared calculation with r || b, solid 
line, and r _L b, dashed line (together with the input distribution, dotted line). The 
calculation with parallel orientation shows a characteristic peak at 6 = |. This corresponds 
to the situation when one end of the dipole, x, is situated in the center of the target where 
the field is strong {U{x.)) ~ and the other end, y, is located in the area where the field 
is very weak U{y) ~ 1. Thus the amplitude A'^ = ^1 — C/"''(x)C/(y)) /N^ ~ 1. On the other 
hand the large dipole perpendicular to the impact parameter axis has both ends in the 
region where the field is not so strong and therefore A^(x, y)_|_ < A^(x, y)|[. 

5.3 Black disc radius and unitarity bound 

The black disc radius Rbd{i~,Y) defines the region in the impact parameter space where 
the amplitude N saturates (for a given dipole size r and rapidity Y). It is defined as the 
solution of the equation 

{N{r,b = RBD,e,Y))g = K, (24) 

with respect to the impact parameter b. As before in Sec. 14.11 in our analysis we choose 
K = 1/2 and average over the angle 6. Thus the black disc radius is a function of the 
dipole size and rapidity. 

The expansion of the black disc area with increasing Y is very important for the 
behaviour of the total dipole-nucleus cross section with energy s ~ sqc^, 

o-(r,y) = 2 / (fhN{v,h,Y). (25) 
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Figure 9: The dipole cross section aBoif-, Y) computed from Eq. (|26]) . plotted as a function 
of rapidity Y for three different dipol sizes r. The dashed hne corresponds to axest = cY. 



This cross section is bounded from below by the cross section integrated over the black 
disc area 

aBD{r,Y) = 2 J (fh N{r,h,Y) e{N - k) ~ 2^i?|^(r,y). (26) 

Thus the problem of the Froissart bound JHl) c(r, 5^) < Y^ for asymptotically high 
energies, can be explicitly studied by looking at the rapidity dependence of the black disc 
radius. 

In Fig. 1^1 we plot ctbd from ()26() as a function of rapidity. For the comparison we 
illustrate the aTest = cY behaviour which would be present if the initial profile in impact 
parameter S{b) ~ exp(— 6^) was preserved. Clearly the behaviour of the black disc cross 
section is much faster than the linear dependence in rapidity. We have found that 

Rbd = Rbd exp(ABD cisY) . (27) 

The extracted value of Abd — 0.6 — 0.7 for r = 0.3 — 1 and A^d — 0.9 — 1 for r = 3 — 30 
at the highest rapidity Y = 11. 
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The value of Xbd increases rapidly for dipole sizes smaller than the given range 
above, mostly due to the pre-asymptotic effects. The value of Xbd is smaller than 
Xbd = 0.87a;jp/(2as) ~ 1.2 which was quoted in Ref. jTUj. The origin of this discrepancy 
can be roughly explained by noting that the power given in |19j most probably overes- 
timates the growth of the black disc radius because it was derived from the analysis of 
the saddle point solution to the linear equation. On the other hand in our numerical 
simulations the value of Xbd = 0.9 — 1.0 for large r is probably closer to the asymptotic 
one. It is due to the fact that for large r, the amplitude has a clear peak for b = r/2 which 
is generated entirely through the evolution. In other words the amplitude in this region 
is less sensitive to the initial profile in b. It is also interesting that this value Xbd con- 
sistent with the relation Abd = 1 /2Xs , based on the following argument of the conformal 
invariance of the equation (see also [MI)- 

Let us take the saddle point solution to the dipole version of the BFKL equation in 
the transverse space jS3 |2S] 

rrn 1 166^ / n rrn 1 \ , , 

N{ro,r,b,Y) ~ -±-—ln--e^p( u^jpY - In^ ^ . (28) 



1662 y3/2 ' \ 1662 14C(3)a,y 

The amplitude is a function of rr^/b"^ only instead of r, 6 separately, which is the conse- 
quence of the conformal invariance. The conditions for the saturation scale (|19|) and black 
disc radius pi]) mean that (roughly) 

. _ 1 Rlo(rX)_^^,^ 



RId(''X) »Q,(>;Y) Q,(b,Y) 

Since in the last equality the right hand side does not depend on the rapidity so should 
not the left hand side too. This means that if we have Rbd ~ exp(ABDas^) and Qg ~ 
exTp{XsasY) then 

Xs = 2Xbd ■ (30) 

Relation ()29() means also that Qs yi, which we have already checked to be approximately 
true (at least for large values of 6 > 7 see Sec. 14. 1|) and also R%b, r. We have verified 
that R%j^{r,Y = 11) ~ with a ~ 0.7 for r = 0.1-1.0 and a ~ 1.0-1.1 for r = 5.0-10.0. 

Clearly the onset of the exponential behaviour in rapidity proportional to exp{XBD(^sY), 
visible in Fig. El signals violation of the Froissart bound. This behaviour is a consequence 
of the power tails in impact parameter N ~ 6^"^, in the same way as the steep exponential 
profile ~ exp(— 62) leads to the linear dependence ~ y of the cross section, compare the 
argument of Heisenberg in j,S8j . The exponential increase of the cross section with rapid- 
ity despite the multiple scattering interactions has been also observed in the Monte Carlo 
study of the amplitude for onium-onium scattering j^B] • 

Since the power-like tail is generated by the long-range contribution, the violation 
of the Froissart bound is caused by the long-range Coulomb-like interactions 19 , which 
in the reality should be suppressed due to confinement. Thus a modification of the BK 
equation by confinement effects is necessary. 
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6 Conclusions 



The solution to the BK equation with the 6-dependence presented in this paper differs 
substantially from the one with translational invariance. The most important result is 
the fact that the exponential profile in b in the initial condition is not preserved by the 
evolution. Instead, the power behaviour is generated for large 6's whose origin comes from 
the structure of the BFKL kernel. Such behaviour leads to the violation of the Froissart 
unitarity bound despite the presence of a local unitarity at fixed impact parameter. The 
violation of this bound is caused by non-suppressed long range contribution, i.e. the lack of 
confinement in the BK equation. This feature of the 6— dependent amplitude is consistent 
with numerical studies of onium-onium scattering with multiple scattering of dipoles pH] . 
and with the qualitative analysis of the BK equation ^U] . 

Thus in order to satisfy the Froissart bound through the evolution one would have to 
modify the evolution kernel in the region of long range contribution in order to incorporate 
confinement. For example, it is evident from Fig. [3 that the naive cut-off - dropping the 
second term in the square brackets in Eq. (|2H|) ~ would do the job. In general, conformal 
symmetry of the integral kernel in the BK equation has to be broken. 

Another interesting feature is the dependence of the solution on the dipole size, which 
shows that the amplitude is saturated in the limited range of scales: l/Qs(b,Y) < r < 
RnibjY). Of course, the true behaviour for large dipole sizes should be modified by the 
long-distance, confinement physics. 
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